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Abstract 

^ I The paper presents derivation and interpretation of one type of 

CN ■ variable order derivative definitions. For mathematical modelling of 

L^ . considering definition the switching and numerical scheme is given. 

ly^ I The paper also introduces a numerical scheme for a variable order 

derivatives based on matrix approach. Using this approach, the iden- 
(^ . tity of the switching scheme and considered definition is derived. The 

^^ I switching scheme can be used as an interpretation of this type of def- 

inition. Paper presents also numerical examples for introduced meth- 
ods. Finally, the idea and results of analog (electrical) realization 
^ ! of the switching fractional order integrator (of orders 0.5 and 1) are 

presented and compared with numerical approach. 



1 Introduction 

Fractional calculus is a generalization of traditional integer order integration 
and differentiation actions onto non-integer order fundamental operator. The 
idea of such a generalization has been mentioned in 1695 by Leibniz and 
L 'Hospital. In the end of 19th century, Liouville and Riemann introduced 
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first definition of fractional derivative. However, only just in late 60' of the 
20th century, this idea drew attention of engineers and mathematicians. The- 
oretical background of fractional calculus can be found in [U |2l |3l HI [5] . Basic 
analysis of fractional differential equations was presented in [HI [71 [8] . Frac- 
tional calculus was found a very useful tool for modelling behavior of many 
materials and systems, especially those based on the diffusion processes. One 
of such devices that can be modeled more efficiently by fractional calculus 
are ultracapacitors. Models of these electronic storage devices, whose capac- 
ity can be even thousands of Farads, based on fractional order models were 
presented in P [lOl [IH [I2] . 

Recently, the case when the order is changing in time, started to be 
intensively developed. The variable fractional order behaviour can be met for 
example in chemistry (when the properties of the system are changing due to 
chemical reactions), electrochemistry, and others areas. In ^3\, experimental 
studies of an electrochemical example of physical fractional variable order 
system are presented. In p3], the variable order equations were used to 
describe a history of drag expression. Papers [121 [13 [IS] present methods for 
numerical realization of fractional variable order integrators or differentiators. 
The fractional variable order calculus also can be used to obtain variable 
order fractional noise [17], and to obtain new control algorithms [18]. Some 
properties of such systems are presented in [19] . In the paper [20] , the variable 
order interpretation of the analog realization of fractional orders integrators, 
realized as domino ladders, was presented. The applications of variable order 
derivatives and integrals can be found also in signal processing [5]. 

The rest of the paper is organized as follows. Section [2] presents existing 
generalizations of Grunwald-Letnikov definition of fractional order deriva- 
tives. In Section [31 a switching scheme for practical implementation of vari- 
able order derivative is given and studied. Section [3] presents also a general- 
ization of the matrix approach for switching order and derivation of identity 
of the switching scheme and the second type of definition. Section [H presents 
numerical examples of the proposed methods compared to the analytical so- 
lutions. Finally, Section [S]presents an analog realization of the switched order 
integrator and comparison of obtained results to the numerical solutions. 



2 Fractional variable order Grunwald-Letnikov 
type derivatives 

As a base of generalization onto variable order derivative the following defi- 
nition is taken into consideration: 

Definition 1. Fractional constant order derivative is defined as follows: 

r=0 ^ ^ 

where n = \t/h\ . 

According to this definition, one obtains: fractional derivatives for a > 0, 
fractional integrals for a < 0, and the original function f{t) for a = 0. For 
the case of order changing with time (variable order case), three types of 
definition can be found in the literature [21], |22]. The first one is obtained 
by replacing of constant order a by variable order a{t). In that approach, all 
coefficients for past samples are obtained for present value of the order and 
is given as follows: 

Definition 2. The 1st type of fractional variable order derivative is defined 
as follows: 



r=0 



In Fig. [U plots of unit step function l(t) derivatives (according to Def. [2]) 
are presented for ai{t) = —1, a2(t) = —2, and 

f-1 forO<t<l, 
"'<" = {-2 forli<<2. « 

The second type of definition assumes that coefficients for past samples 
are obtained for order that was present for these samples. In this case, the 
definition has the following form: 

Definition 3. The 2nd type of fractional variable order derivative is defined 
as follows: 

r=0 ^ ^ 




Figure 1. Plots of unit step function derivatives with respect to the 1st type derivative 

(given by Def. ^ 

In Fig. El plots of unit step function derivatives (according to Def. [3]) 
are presented for ai(t) = —1, 0^2 (^) = —2, and 0:3 (t) given by ([1]). The 




Figure 2. Plots of unit step function derivatives with respect to the 2nd type derivative 

(given by Def. ^ 

third definition is less intuitive and assumes that coefficients for the newest 
samples are obtained respectively for the oldest orders. For such a case, the 
following definition applies: 

Definition 4. The 3rd type of fractional variable order derivative is defined 
as: 

{-ly fa{rhy 






r=0 



fit-rh). 



In Fig. [31 plots of unit step function derivatives (according to Def. H]) are 
presented for ai(t) = —1, a2{t) = —2, and a^it) given by ([I]). 




Figure 3. Plot of unit step function derivatives with respect to the 3rd type derivative 

(given by Def . S]) 



From the comparison of the three types of derivatives, it can be seen 
from the plots presented above that there are crucial differences between 
the nature of derivatives of switching order (for constant order derivatives 
all the definitions yield the same behavior). Namely, in the case of the 1st 
type derivative, at the switching time instant, the derivative output "jumps" 
from the plot of constant order ai to the plot of constant order 02 and keeps 
to follow it. In the case of the 2nd type derivative, at the switching time 
instant, the derivative output stops fitting the plot of constant order ai and 
starts to integrate like at the beginning of the plot of constant order a2 but 
starting from another initial value. Concerning the 3rd type derivative, at the 
switching time instant, the plot of switching order derivative stops fitting the 
plot of constant order ai and does continue to integrate with constant order 
a2 (without "jumping" to the plot of 0:2 - unlike the 1st type derivative) like 
it is for the plot of constant order 02 at the same time interval, i.e., for t > 1. 



3 Practical implementation and numerical scheme 
of switched (variable) order derivative 

In this section, the routines and schemes for switching order derivative are 
presented. For simplicity, we start with the simplest case of order switching, 
namely switching between two real arbitrary constant orders, e.g., ai and a2- 
Next, this idea will be generalized for a multiple-switching (variable order) 
case. 



3.1 Simple-switching order case 

The idea is depicted in Fig. HI where all the switches Si, i = 1,2, change 
their positions depending on an actual value of a{t). If we want to switch 
from ai to 02, then, before switching time T, we have: Si = b, S2 = a, and 
after this time: Si = a and S2 = b. At the instant time T, the derivative 
block of complementary order 0:2 is pre-connected on the front of the current 
derivative block of order ai, where 



Ci2 = 0(2 — Cil- 



(2) 



If a2 < 0, then tD"^ corresponds to integration of f(t); and, if 02 > 0, then 
tD^^ corresponds to derivative of f{t), with appropriate order a2- 
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Figure 4. Structure of simple-switching order derivative (switching from ai to 02) 



Now, the numerical scheme corresponding to the above derivative switch- 
ing structure is introduced. The matrix form of the fractional order derivative 
is given as follows [231 [2i 
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number of samples. 



— i^a , and h = t/k is a time step, A; is a 



Lemma 1. For a switching order case, when the switch from order ai to 
order 0:2 occurs at time T, the numerical scheme has the following form: 
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The order 02, appearing above, is given by relation ([2]). 

Proof. The signal incoming to the block of derivative ai can be described as 
follows 

/ /(O) \ / /(O) \ 
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Until time T, the input of ai-block obtains the original function /(t), so in 
matrix VF(q;2, k., T) we have an identity matrix. From time step T the input 
signal is passes through the block of derivative 0:2 and we have the sub-matrix 
W{a2, k — T) that is responsible for starting the 0.2 derivative action from 



time T. That signal passes to the ai-block and has the following matrix 
form: 
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Example 1. Consider integration of the unit step function f(t) = l(t), with 
the switching variable order a{t) taking value ai = —luntil switching time 
T = 4, and after this time, order a2 = —2. 

Numerical calculations performed for k = 5 yield: 
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where a^ 



-1. 



3.2 Multiple-switching (variable order) case 

In general case, when there are many switchings between arbitrary orders, we 
have the following structure, presented in Fig. O When we switch from the 
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Figure 5. Structure of multiple-switching order derivatives 



order cxj^i to the order Uj at the switch-time instant Tj_i, for j = 2, 3, 
we have to set: 

a for i = 1, . . . , j — 1, 

b for i = j, 

and the pre-connected derivative block (on the front of the previous term) is 
of the following complementary order: 



5,: 



where 
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The numerical scheme describing the already presented general case of struc- 
ture allowing to switch between an arbitrary number of orders is given be- 
low. What is very important, the numerical scheme of multiple-switching 
case (when order is switched in each sample time) is equivalent to the 2nd 
type of variable order derivative. 

Theorem 1. Switching order scheme presented in Fig. is equivalent to the 
2nd type of variable order derivative (given by Def. [^. 

Proof. Let us assume that the order of derivative changes with every time 
step, which gives a variable order derivative, and is given as follows: 

k 

ak = J2^J^ (4) 



3=0 



where, in this case, ao = ao is a value of initial order. Using Lemma [T] the 
following numerical scheme is obtained 
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The first switching matrices can be described as the following block ma- 
trices: 
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is a vector with coefficients given by ^. For a switching in the next sample 
time, we obtain the following numerical scheme: 
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In the case oi k — th switchings of order, we have the following form of the 
switching matrix W{ao, k, 0)W{di, k,l) ■ ■ ■ W{ak, k, k), i.e.. 
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or shortly, using the relationship 
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The coefficients of the matrix above are identical with the coefficients given 
by Def. [3l which completes the proof. D 
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4 Numerical examples 

This section contains numerical examples of the proposed methods, com- 
puted in Matlab/Simulink environment, using the dedicated numerical rou- 
tines [25], developed by the authors. 

Example 2. Let us consider a variable order integration of unit step function 
fit) = l(t) with the following sequence of integer orders A = {— 1, — 2, — 3, — 1} 
switched every one second, i.e., 



a{t) 



' -\ for 0<t<l, 

-2 for l<t < 2, 

-3 for 2<t < 3, 

-1 for 3<t< 4. 



The analytical solution of such integration is the following 
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dr 

for < t < 1, 

for 1 < t < 2, 

for 2 < t < 3, 

for 3 < t < 4. 
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Figure 6. Plot of the variable order integral of unit step function (Ex. [2]) 
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Figure 7. Errors between analytical and numerical solution for different integration 

steps (Ex. [2) 



The plot of integration result is presented in Fig. O The difference between 
analytical and numerical solution, for different integration steps h, is depicted 
in Fig. [^ 

Example 3. Let us consider a variable order integration of constant function 

f{t) = l(t) with the following sequence of fractional orders A = {—0.4, —1.8, —0.5, —2.5} 

switched with every one second, i.e.. 



a{t) 



'-0.4 for < t < 1, 

-1.8 for 1 < t < 2, 

-0.5 for 2 < t < 3, 

1^-2.5 for 3 < t <4. 
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The analytical solution of such integration is the following 

' Di for < t < 1, 

oDfV(t)-r' /orl<t<2, 
I ^3 /or 2 < t < 3, 
^Di for 3 < t < 4, 



where 



Di = 1.127t°•^ 

D2 = 0.597 {t-iy-^ 
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D3 = l.USVt^^ + 0.597 ((t - 1)^-^ - (t - 2)^-^) 

+ 1.127 (t°-^-(t-l)°-^), 
D4 = 0.3(t - 3)'-^ + 0.597 ((t - 1)^-^ - (t - 2)^-«) 

+ 1.128 {Vt^ - Vt^) + 1.127 (t°-' - (t - 1)°-') . 
r/ie plot of integration result is presented in Fig. 0. 
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Figure 8. Plot of the variable order integral of constant function (Ex. [3]) 



5 Analog realization of the second type of 
fractional variable order derivative 

5.1 Analog realization of the half order integral 

In this paper, the following method of half order integrator implementation, 
introduced in [26] and meticulously investigated in [201 i2Z] , will be used. The 
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scheme of this method is presented in Fig. [HI 

i{t) -^1 -^2 Ri 



u{t) 



Ci: 



^2 



Ci: 



Ci: 



Figure 9. Domino ladder type analog model of half order integrator 

Based on the observation, the values of Ri and Ci are chosen in order 
to satisfy the required low frequency limit. Value i?2 is chosen in order to 
satisfy the condition Ri/R2 ^ 0.25, and the realization length m (number 
of branches) is chosen in order to satisfy required frequency range of the 
approximation. 

5.2 Experimental setup 

An analog realization of switching system directly based on the switching 
scheme given in Fig. H] is presented in Fig. [TOl This scheme can be easily 
extended to the switching rule presented in Fig. |5l The experimental setup 
was prepared in order to modelling simple-switching case. The structure of 
the half order integrator, used in the setup, is based on the domino ladder 
approximation presented in Section 15.11 However, the first order integrator 
is realized according to traditional scheme based on capacitor. Because, 
in both cases of realizations with operational amplifiers, the signals with 
inverted polarization were obtained. This required amplifiers with gain —1 
for each integrator (amplifiers A2 and A4). Integrators based on amplifiers Ai 
and ^3 contain gain resistors R = 9.65kQ and impedances Zi and Z2 chosen 
according to the type of realized orders. As a realization of switches 5*1 and 5*2 
integrated analog switches DG303 were used. In order to obtain impedance 
order equal to —0.5, the domino ladder approximation with the following 
elements Ri = 2kQ, R2 = 8.2kQ, Ci = 470nF and of implementation length 
771 = 100, was used. The experimental circuit was connected to the dSPACE 
DS1104 PPG card with a PC. 
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Figure 10. Analog realization of the 2nd type of fractional variable order integral 



5.3 Experimental results 

5.3.1 Switching between order a = —0.5 and a 



-1 



For a case of switching between order a = —0.5 and a = — 1, in the structure 
presented in Fig. [T0|, both impedance Zi and Z2 contains domino ladders 
structures each of order —0.5. This implies that operational amplifiers Ai 
and ^3 realize half order integrators circuits. Before switching, only one half 
order integrator is used and the whole system possesses —0.5 order. After 
switching, both half order integrators are connected in series, which results 
that the system order is equal to —1. 

The experimental results of integrator with switched orders from a = 
—0.5 to a = — 1 compared with numerical results are presented in Fig. [TTl 
The difference between numerical realization of the 2nd type derivative def- 
inition and its analog implementation is presented in Fig. [121 The sam- 
ple time for all measurements was chosen as 0.01 sec and input signal was 
u{t) = 0.01 ■ l(t). The parameters of analog models were obtained by identifi- 
cation based on time domain responses, separately for both orders. Obtained 
transfer function for the case of order a = —0.5 (single half order integrator) 
is 
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and for the case of order a = —1 (both half order integrators in series con- 
nection) 

G2(S) ^ 



0.041s' 



The switching time was equal to 1 sec. 




analog model 

numerical Implementation of 2nd definition 



Figure 11. Results of analog and numerical implementation of switching order 
derivative from a = — 0.5toQ: = — 1 




Figure 12. Difference between analog and numerical implementation of switching order 

derivative from q; = — 0.5toa = — 1 



5.3.2 Switching between order a = —0.5 and a = —1.5 

In this case, impedance Zi is the first order capacitor, and impedance Z2 is a 
—0.5 order domino ladder. The experimental results of integrator with order 
switching from a = —0.5 to a = —1.5, compared to the numerical results. 
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are presented in Fig. [121 The difference between numerical realization of the 
2nd type derivative definition and its analog implementation is presented in 
Fig. [HI By identification the following models were obtained. For order 
a = -0.5 

GAs) ^ 



and for order a = —1.5 



G.(s) 



0.68s0-5^ 
1 



0.0295si-5' 
The switching time was equal to 1 sec. 
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numerical implementation of 2nd definition| 




Figure 13. Results of analog and numerical implementation of switching order 
derivative from a = —0.5 to a = —1.5 




Figure 14. Difference between analog and numerical implementation of switching order 

derivative from a — —0.5 to a = —1.5 



17 



5.3.3 Switching between order a = —1 and a = —1.5 



In this case, impedance Zi is a domino ladder of order —0.5, and impedance 
Z2 is the first order capacitor. The experimental results of integrator with 
order switching from a = —1 to a = —1.5, compared to the numerical results, 
are presented in Fig. [T3 The difference between numerical realization of the 
2nd type derivative definition and its analog implementation is presented in 
Fig. [161 By identification the following models were obtained. For order 
a = —1 



and for order a 



■1.5 



G.is) 



0.038s' 



1 



0.027si-5 
The switching time was equal to 0.7 sec. 
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—numerical Implementation ol 2nd definition 



Figure 15. Results of analog and numerical implementation of switching order 
derivative from a = — Itoa^— 1.5 



6 Conclusions 



The paper presented a switching scheme for the second type of variable or- 
der derivative. The numerical scheme, based on matrix approach, for that 
switching scheme was introduced and investigated. It was shown that ob- 
tained numerical scheme is equivalent to the 2nd type of fractional variable 
order derivative. This switching scheme can also be used as an interpre- 
tation of the second type of the definition. It is also worth to notice that 




Figure 16. Difference between analog and numerical implementation of switching order 

derivative from a=:— ltoQ! = — 1.5 



the rule given by the switching scheme can itself be a definition of a vari- 
able order derivative and the relation given by Def. |3]is just a consequence 
of this scheme. Presented interpretation allows to better understand be- 
haviour of the definition and, in general, switching process in variable order 
systems. Based on this, it can give rise to more appropriate choice of def- 
initions type, depending on particular application. Additionally, an analog 
modelling method of switched order integrator was introduced and exam- 
ined. Obtained results were compared with the numerical implementation 
and show high accuracy of the introduced method. 
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